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The percolating phase of whatever random percolation process resembles the confining vacuum of a 
' gauge theory in most respects, with a string tension having a well-behaved continuum limit, a non 

trivial glueball spectrum and a deconfinement transition at a well determined temperature T c . Simple 
| numerical experiments reveal an underlying, strongly fluctuating, confining string, with an internal 

, vortex structure formed by a core trapping inside a Coulomb-like phase composed by the vacuum at 

the percolation threshold. The width of the core almost coincides with 1 /T c and it turns out to be 
£^ ■ separated form the confining vacuum by a domain wall of definite thickness. 

a 

1 Introduction 

Recently, it has been pointed out that three-dimensional percolation can be thought as a 

m 

gauge theory in disguise . Although percolation is a purely geometrical process in which 
no dynamics is involved, it captures all the salient features of a confining vacuum. 

Take for instance a bond percolation process. The ensemble of configurations is ob- 
tained simply by populating each of the links of a 3D lattice independently with occu- 
pation probability p. The connected graphs made by occupied links are called clusters. 
When p exceeds a threshold value p c there is a percolating cluster in the infinite lattice. 

The key ingredient to extract from this ensemble the relevant information on the un- 
derlying gauge theory is the definition of the percolation counterpart of the Wilson oper- 
ator Wj, that we associate to whatever loop 7 of the dual lattice. The rule of the game is 
very simple. We set W 7 = 1 if there is no cluster of the configuration topologically linked 
-C ■ to 7, otherwise we put W 1 = 0. Therefore, the vacuum expectation value (Wy) coincides 

with the average probability that there is no path in any cluster linked to 7. As in usual 
gauge theories, evaluating these quantities yields the main physical properties of this sys- 
tem. In this way it has been shown that in the percolating phase the theory i s conf ining, 

in analogy to what happens in gauge theories with percolating center vorticeP^^! The 
string tension a and the other physical observables have the expected scaling behaviour 
dictated by the universality class of 3D percolation, therefore such a theory has a well 

defined continuum limit. Moreover it has a non-trivial glueball spectrum^! and a second- 
order deconfining transition at finite temperature T c with a ratio T c j ' sfa ~ 1.5 which 
turns out to be universal, i.e. it does not depend on the kind of lattice utilised nor on the 

specific percolation process considered ( site or bond percolation)^. 

A particularly interesting feature is the occurrence of universal shape effects in the 

Wilson loopu^l which are usually ascribed to the quantum fluctuations of a string. Then 
the question arises: where and how are encoded the relevant degrees of freedom of the 
underlying confining string? This is the main issue discussed in this paper. The answer 
we find, though surprising, confirms once again that percolation is a full-fledged gauge 
theory: the confining string reveals a non-trivial microscopic structure which is believed 
to be shared by whatever confining gauge theory, according to the picture of the confining 

vacuum as a dual superconductor^^. 



Figure 1 : Distribution of the average occupancy probability in the plane of a rectangular 10 x 20 Wilson 

loop at p — 0.26. 



The great advantage of percolation is that its simplicity allows to explore regions that 
are still inaccessible to the other gauge systems from a computational point of view. 

It turns out that the confining string in percolation is made of a core in a Coulomb- 
like vacuum separated from the confining vacuum by a domain wall. The width 1Z of 
the core coincides almost exactly with the inverse of the deconfining temperature, while 
the thickness of the domain wall is about TZ/2. This confining string is not stuck to the 
minimal surface bounded by 7, but strongly fluctuates, sweeping a large volume which 
grows with the size of the Wilson loop. 

2 Wilson loop as a source of the flux tube 

Each bond configuration of a percolating process in a cubic lattice can be mapped into 
a plaquette configuration of the dual lattice by setting empty any plaquette orthogonal to 
an occupied dual link and vice versa. This mapping is one to one, therefore the question 
whether there is a closed path in a cluster that is linked with a loop 7 of the dual lattice 
is equivalent to the existence of a surface S (at least) of occupied plaquettes having 7 
as boundary, i.e. = 7. The cluster of this kind of surfaces can be considered as 
the outgoing flux of the "gauge" field generated by a point-like source describing the 
world-line 7. 

The crucial observation which allows to extract a lot of information form the apparent 
white noise of the percolation ensemble is that not all the configurations are compatible 
with the existence of such a source, but only those with W 1 = 1. In this subset the 
distribution of occupied plaquettes is not at all flat and may be used to study the spatial 
distribution of the gauge field. When p is below the percolating threshold the flux is not 
very constrained and is spread out in the whole space. On the contrary, in presence of a 
percolating cluster the volume at the disposal of the gauge flux is much more tight, thus 
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Figure 2: Profile of the surface drawn in the previous figure along the line at y = 10. The horizontal straight 
line denotes the threshold value of bond percolation. 



the percolating cluster acts as a sort of superconducting medium, squeezing the gauge 
flux in a thin tube. 

It is useful to associate to each link t of the lattice a projector n(£) which is 1 when £ 
is occupied an otherwise. We have trivially (tt(£)} e p = P, where E p is the whole en- 
semble of the percolation process. Much less trivial is the average occupancy probability 
of the links in the subset of configurations with W 1 = 1, that we denote as 
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(7T(£)) {Wj = l } , (D 



and we use it as a local probe to explore the spatial distribution of the flux. As an example, 
the quantity -kJ as a function of the links £ orthogonal to the plane of a 10 x 20 rectangular 
Wilson loop in the percolating phase at p = 0.260 is drawn in Fig{2 The data are taken 

using a simple burning epidemic type of algorithm^. The sinking along the perimeter of 
the rectangle is an effect caused by the finite clusters: since the clusters near 7 have a 
finite probability to be linked to it, the configurations selected by W 1 = 1 have a cluster 
density smaller with respect the whole ensemble E p a . Therefore an observer near 7 
experiences an occupation probability 7r 7 < p. One has even 7r 7 < p c , as Fig|2] shows. 
Hence the vacuum around the 7 source is in a non-percolating phase. This important 
fact will be the starting point of the next section. 

Inside the rectangular surface encircled by 7 the average probability it] reaches a 
plateau p < p because of the presence of the infinite, percolating cluster. Actually this 
very condition is sufficient to assure that the vacuum in which the Wilson loop is embed- 
ded is confining. Indeed it is easy to prove the following exact inequality for the string 
tension 
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"Also the deepening at the corners of the rectangle is a finite cluster effect, because the linking probability 
of a cluster is proportional to its overlap with 7 and such an overlap is larger near a corner. 
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Figure 3: A sketch of the microscopic structure of the flux tube. The Coulomb phase inside is made of 
vacuum at the percolation threshold p c , while the outside vacuum is at p > p c . 



To prove it, note that the probability of finding k occupied links among the N links 
orthogonal to a rectangle 7 = r x t is binomial by construction b : : 
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In our case N = rt and we are assuming that k ~ pN with p < p. For large k and N 
we have 
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On the other hand, not all the configurations with p ~ are compatible with W r ,t = 1> 
of course, hence we get the exact inequality 



(W r ,t) <Pn(PN) , N = rt, 



(5) 



which yields, in turn, eq.©. 



3 The intrinsic width of the flux tube 

What exactly is the confining string made of? The distribution of the occupied links in 
close proximity of the source 7 of the flux tube is much smaller than its asymptotic value 
and even smaller than the percolation threshold, as Fig. \\\ and El clearly show. This sup- 
port the idea that the confining string may have a microscopic structure similar to that of a 
(dual) Abrikosov or Nielsen-Olesen vortex , with a core made of a Coulomb-like vacuum 
as contrasted with the surrounding medium, which is in the confining phase. This core 
conveys the whole flux in a tube of small but non-zero thickness 1Z C (see Figl^J). Some 
nume rical an d analytical tests of this scenario were performed long ago on the Z2 gauge 

modeinilZI The semi-classical descriptions of these extended objects neglect an impor- 
tant quantum property of the confining string. In fact, in order to account for the universal 
shape effects of large Wilson loops (W 7 ) one is forced to assume that the string world- 
sheet belongs to a rough phase, hence it should resemble much more closely a strongly 

''With a similar reasoning one could "prove" quark confinement out of percolation of center vortice/. 
There is however a flaw in the argument, because center vortices, contrary to the links of random percolation, 
are correlated, therefore they cannot obey a binomial distribution. 

c For an updated discussion on this argument and a complete list of references seel^l 



fluctuating fluid interface rather than a static, smooth structure. As a consequence, we do 
not expect to find a flux tube stuck to the minimal surface encircled by 7. Nonetheless, 
the average occupancy probability in the whole space would provide clues to its internal 
structure. In this way, enlarging the size of the lattice, one discovers that the confining 
string sweeps a very large volume surrounding the loop 7. More precisely, we found that 
the sum 

over the n links I of the lattice A is less than the naively expected value p n and their 
difference pn — w approaches a constant value depending only on the size of 7. This can 
be considered as a first hint on the finiteness of the world-volume of flux tube. We start 
by considering a rather crude approximation of the flux tube, by partitioning the n links 
of the lattice into two subsets 

n = m + n e , (7) 

where rij is the number of links lying inside the Coulomb-like core of the string and n e 
the number of those lying outside, in the confining vacuum. 

The average number of occupied links in the confining phase is pn e , while the links 
belonging to the core are occupied with a smaller probability p < p c in order to prevent 
linking. This fact is manifest near the fixed boundary 7, where the string cannot fluctuate 
too much, as shown in Fig|2] Thus we can rewrite eq.© as 

w =prii+pn e , (8) 

which tells us that the volume of the core can be expressed simply as 

n t = P -^. (9) 
p — p 

In order to estimate n; L we need to know the value of p. A simple scaling argument sug- 
gests p = p c . When p — > p c , the flux tube is expected to grow because the infinite, 
percolating cluster is crumbling, then the flux can spread out in the whole space. There- 
fore, from one hand we expect in this limit ni ~ n. On the other hand in this very limit 
the confining string disappears, hence we should have w ~ pn. According to eq.Q, 
these two conditions are consistent only if p = p c . 

Assuming the world- volume of the flux tube to be a slab of size A TZ C , with A the area 
of the minimal surface encircled by 7, we are led to the following effective definition of 
the intrinsic thickness of the flux tube 

A A(p - p c ) 

Some numerical estimates of 1Z C are reported in Tab^ They seem to approach asymp- 
totically an approximate scaling behaviour, as expected for a physical quantity. Though 
it gives a length scale of the order of the transverse size of the flux tube, it should not be 
taken too literally. First, the splitting of the links in two subsets of those inside or outside 
the Coulombic core is an over-simplification: there should be a domain wall between 
the two vacua with an average occupancy probability that interpolates between p and p c - 
This would produce an enlargement of the flux tube which is not taken into account in 
eq.fllOli. Similarly, the number of interior links does not coincide exactly with the volume 



Figure 4: The slab geometry to probe the internal structure of the flux tube. The transverse direction is 
periodic and of size L±. 7 denotes a rectangular Wilson loop orthogonal to the compactified dimension. 
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Figure 5: The intrinsic width of the confining string as a function of the transverse dimension L± of the 
periodic slab in which is trapped. The data are taken at p = 0.26 with a spatial square Wilson loop of size 
£ = 15. The horizontal arrow denotes the asymptotic value R c (oo) of the intrinsic width. 



Table 1: The intrinsic width 1Z C of the flux tube generated by a square Wilson loop of side L. The data 
are taken at two different values of p. The corresponding values of a and of the inverse deconfinement 
temperature 1/T C are also reported. All the physical quantities are expressed in lattice spacing units. 
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Figure 6: Vacuum expectation value of a square Wilson loop of size £ — 15 at p — 0.26 as a function of the 
transverse periodic dimension L± of the lattice. The arrow indicates its asymptotic value. 



of the core, the basis of the slab occupied by the world- volume of the flux tube is actually 
larger than A, and so on. Nevertheless we can still extract further information on the 
microscopic structure of the confining string combining these results with another kind 
of numerical experiments. 

Lo ng tim e ago it was observed that in three-dimensional S'[/(2)Eland Z2 gauge 

models^ ^ the intrinsic width of the flux tube almost coincides with the inverse of 
the deconfinement temperature (see also 1 ^. This result was obtained by studying the 
response of the vacuum expectation value of a large, planar Wilson loop when one varies 
the size L± of the dimension orthogonal to the loop in a periodic lattice, like in Fig|4] 

We repeated such a numerical experiment in percolation, using the same setting and 
measuring the intrinsic radius 1Z C as defined in eq. dTOl . The surprising results are reported 
in Fig|5]and|6] Let us follow the behaviour of R C (L±) starting from L± large enough, 
where R c coincides with its asymptotic value 1Z C = R c (oo). As L± decreases, we 



observe first an apparent growth of R C (L±) starting at a value of L± slightly larger than 
the asymptotic value 1Z C , the reason being perhaps the underestimated contribution of the 
domain wall. The apparent increasing of R C (L±) should be ascribed to the fact that the 
two domain walls separating the core from the confining vacuum are pushed one against 
the other squeezing a fraction of the core. When L± is of the order of 1/T C , the external 
confining vacuum disappears and the domain walls melt in the Coulomb-like core. At 
this point the flux tube fills the entire interval L± and loses its own width, which now 
coincides, for L± < l/T c , with the transverse size of the lattice. This fact is corroborated 
by value of the slope in this region, which is about 1 . 

In conclusion, the response of the system to the variation of the size of the trans- 
verse dimension seems to confirm the scenario of a confining string made of a core in 
a Coulomb-like vacuum separated from the confining medium by a domain wall. The 
width 1Z of the core almost coincides with the inverse of the deconfining temperature, 
like in other gauge systems, while the thickness of the domain wall seems to be about 
K/2. 
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